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Abstract 

Recent "work on the behaviour of locahsed states in pattern forming partial differential equa- 
tions has focused on the traditional model S"wift-Hohenberg equation "which, as a result of its 
simplicity, has additional structure — it is variational in time and conservative in space. In 
this paper we investigate an extended S"wift-Hohenberg equation in "which non-variational and 
non-conservative effects play a key role. Our "work concentrates on aspects of this much more 
complicated problem. Firstly wc carry out the normal form analysis of the initial pattern form- 
ing instability that leads to small-amplitude localised states. Next "wc examine the bifurcation 
structure of the large-amplitude localised states. Finally we investigate the temporal stabil- 
ity of one-peak localised states. Throughout, we compare the localised states in the extended 
S"wift-Hohenberg equation "with the analogous solutions to the usual S"wift-Hohenberg equation. 

Key"words: homoclinic snaking, pattern formation, dissipative solitions 



1 Introduction and background 

The study of the spontaneous emergence of patterns of activity out of homogeneous states has a 
long history, motivated by a "wealth of examples in fluid and solid mechanics, and more recently 
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extended to nonlinear optics, granular media, chemical reactions and mathematical biology. General 
surveys are given, from different viewpoints, by Cross and Hohenberg [10], Cross and Greenside [9], 
Hoyle [16] and Pismen |23j . 

In many cases the stable patterns that appear are roughly periodic in space and extend through- 
out the bulk of the experimental or computational domain, growing smoothly in amplitude as a 
control parameter increases above a critical value. In other cases the system displays hysteresis: the 
pattern appears abruptly at finite amplitude as the system undergoes bifurcation, and the pattern 
persists as the control parameter is then reduced below the bifurcation point. This latter case is 
often referred to as the 'subcritical' case, in contrast with the former which is the 'forward' or 
'supercritical' case. 

One of the most popular model equations for the examination of the dynamics of pattern- 
forming systems of this kind on a domain C M is the Swift-Hohenberg equation [25] in one 
spatial dimension with quadratic-cubic nonlinearity: 



where u{x, t) is a scalar variable which describes the pattern forming activity, and r and b are real 
control parameters. Typical analyses [3] fix 6 and treat r as the primary bifurcation parameter. The 
trivial state u{x, t) = is linearly stable in r < and undergoes a pattern forming instability at r = 
0, where the spatial coordinate in ([TJ is scaled so that the initial instability is to perturbations with 
wavenumber k = 1. In r > the trivial state is unstable to steady, spatially periodic perturbations 
with a range of wavenumbers surrounding k = 1. The secondary bifurcation parameter b determines 
the criticality of the pattern forming instability at r = 0: it is supercritical if 6^ < 27/38 and 
subcritical if 6^ > 27/38. 

Equation ([T]) has several important symmetries. Firstly, it is equivariant under spatial reflec- 
tions {x,u) — )• {—x,u). Secondly, it is equivariant under the following inversion involving the 
parameters: {x, u; r, b) — )• (x, —u; r, —b). In consequence, the behaviour of ([1]) can be fully classified 
by considering only the 6 > half of the parameter plane. 

We note that a common variant of ([1]) replaces the quadratic-cubic nonlinearities bu^ — with 



dtu = ru — {1 + dxx)^ u + bu^ — v? 
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the cubic-quintic combination cu^ — ti^ [6]. This results in a subcritical instabihty for all c > 0, 
and is appropriate if there is an additional symmetry u — )• —u of the system, as for example in the 
case of Boussinesq thermal convection with identical upper and lower boundary conditions [16] , 
Of course ([T]) also includes this extra symmetry at 6 = 0, and this is a popular choice to model a 
supercritical pattern forming system. In this article, we are interested in the more generic case in 
which the additional u ^ —u symmetry is absent. 

The dynamics of ([T]) are strongly influenced by the fact that it is variational — i.e., it can be 
written in the form 

dtu = ^ , (2) 

ou 

where the Lyapunov functional J-\vi\ (which we refer to as a free energy) is given by: 

on the bounded domain with suitable boundary conditions on dO. (for example, Neumann). It 
follows that 

^^J^[u] = - l^{dtufdx<0, (4) 

so that the free energy decreases in time along trajectories. Straightforward calculations show that 
J-[u\ is in addition bounded from below. Hence the variational property guarantees that solutions 
converge to equilibria: sustained oscillations, travelling waves, and temporal chaos can not arise. 
The resulting equilibrium profiles u{x) satisfy the time-independent version of ([1]). An equivalent 
description of these equilibria is as trajectories in x of a fourth-order spatial dynamical system. 
For example, a stationary, uniform amplitude pattern of wavenumber k corresponds to a periodic 
orbit in x with period 27r//c. The spatial dynamical system thus derived is conservative, and the 
quantity 

H = -^ir- l)n2 + (d^uf - ^{d^^uf + (a,u)(9,,,u) - hu' + ^u^ (5) 

is independent of x. Both of these properties — variational in t and conservative in x — aide in 
the analysis of ([T]). 

On spatially extended domains = M, the Swift-Hohenberg equation ([1]) exhibits a variety 
of equilibria beyond the uniform amplitude patterns mentioned above. In the subcritical regime, 
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this includes stationary, spatially localised states for which the amplitude of the pattern decays 
to u —7- as X —7- iboo. In the spatial dynamical system, these localised profiles correspond to 
orbits homoclinic in x to the trivial fixed point n = 0. Recent short reviews of the generation 
of localised states in this context are given by Knobloch [T^] and Dawes The bifurcation 

analysis of the initial pattern-forming instability at r = (equivalently referred to in the context 
of fourth-order dynamical systems as a reversible 1:1 resonance or a Hamiltonian-Hopf bifurcation) 
was initiated by looss and Peroueme [15] and developed by later authors, notably Woods and 
Champneys |26J, Coullet et al. [8], and Burke and Knobloch [5]. The normal form analysis shows 
that in the subcritical regime a family of steady, small amplitude localised states of the form 



bifurcates from n = into r < 0, along with the uniform amplitude patterns. Terms present 
in the Swift-Hohenberg equation that are manifested beyond all algebraic orders in the normal 
form break the normal form symmetry associated with cj) and select (p = and = vr as the 
only physical solutions to ([T]) — see Refs. [71 [20]. The profiles along these two branches are even- 
symmetric under spatial reflection (i.e., they are invariant under the reversibility transformation). 
The branch associated with (p = includes profiles with a local maximum in u at the midpoint, 
and the branch associated with (p = tt includes profiles with a local minimum in u at the midpoint. 
We refer to the former branch as LO and the latter as LI. These two branches of even-symmetric 
localised states persist to finite amplitude where they undergo homoclinic snaking — a sequence 
of saddle-node bifurcations that cause the branches to intertwine as they oscillate back and forth 
across a parameter range called the snaking or pinning region [H [26] . The resulting bifurcation 
structure, shown in figure [H also includes a sequence of so-called rung branches which emerge from 
LO and LI in pitchfork bifurcations located near the saddle- node bifurcations and cross-link the 
two snaking branches. The profiles that make up LO and LI resemble several wavelengths of a 
uniform amplitude pattern connected to the trivial state by a symmetrically related pair of fronts. 
The profiles on the rung branches are similar to those on the snaking branches, but they are not 
symmetric. These profiles spontaneously break the reversibility symmetry of ([1]), so each point on 
the rung branches in figure [1] actually represents two different profiles with identical norm which 




(6) 
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Figure 1: (a) Bifurcation diagram of stationary solutions to equation ([T]) at & = 1.8, plotted in terms of 
the norm 11^11^2 = {jQ^'^i^)d-xY^'^ ■ Shading indicates the snaking region. The snaking branches LO and LI 
include even-symmetric localised states; the arrows indicate that on 17 = K the snaking continues indefinitely. 
The rung branches which cross-link the snaking branches are also shown. The branch P of spatially periodic 
patterns satisfies H = 0, and includes the Maxwell point M at which 7^ = 0. The norm of solutions on 
P is rescaled so that this branch can be displayed on the same scale as the branches of localised states. 
Solid/dashed curves indicate stable/unstable solutions, (b) Profiles from several saddle-node bifurcations of 
the snaking branches — profiles (i), (iii) and (v) are from LO, and profiles (ii), (iv) and (vi) are from LI. 

are related to each other by reflection in x. 

The linear stability of the various stationary profiles is determined by linearizing ([1]) about 
the state. The even-symmetric localised states from LO and LI are unstable at onset and change 
stability at each saddle-node bifurcation; profiles from the segments of the snaking branches that 
slant 'up and to the right' on the bifurcation digram in figured] are stable and those that slant 'up 
and to the left' are unstable. All of the asymmetric profiles from the rungs are unstable. 

The variational and conservative properties of the Swift-Hohenberg equation help considerably 
in understanding these localised states and the associated snaking bifurcation structure [6]. For 
example, the fact that the spatial dynamics associated with ([T]) is conservative determines the 
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wavelength (i.e., the spatial period) of the pattern within the localised states. At fixed r, there 
typically exists an entire family of stationary, spatially periodic patterns up{x;k) parameterised 
by the wavenumber k. The particular pattern that is selected to appear within the localised state 
must lie in the level set H = 0. Figure [T] includes the branch P of spatially periodic states defined 
hy H = 0, and the pattern wavenumber k varies with r along this branch to satisfy the H = 
constraint. Careful measurement of the numerically computed localised states confirms that this 
branch of patterns correctly predicts the wavenumber variation k{r) within the localised states, at 
least when the localised states are sufficiently wide. The variational property of ([T|) is also useful 
in understanding the localised states. The free energy of the uniform amplitude patterns varies 
along P. The so-called Maxwell point M is the r value at which the pattern on P has the same 
free energy as the trivial state — i.e., it is defined by the two constraints H = = 0. Away from 
the Maxwell point, there is a free energy mismatch between the pattern within the localised state 
and the trivial background state, with the pattern energetically favoured above the Maxwell point 
and the trivial flat state favoured below. This free energy mismatch is countered by the energy 
associated with the fronts, which pins the fronts to the interior pattern. The stationary localised 
states within the snaking region correspond to critical points of the free energy landscape (local 
minima for the stable localised states, local maxima and saddles for the unstable localised states 
from the snaking branches and the rungs). Outside the snaking region the free energy mismatch 
is sufficiently large to eliminate the critical points, de-pinning the fronts and forcing them to drift. 
The temporal dynamics of localised initial conditions at values of r outside the snaking region 
support this description. The Maxwell point therefore serves as an organizing center for the snaking 
structure and the localised states. 

It is clear that the special properties of the Swift-Hohenberg equation confer additional prop- 
erties which, while extremely useful in understanding the bifurcation structure of localised states, 
may not be valid generically. The generic situation in which Turing instability occurs in one di- 
mension is spatially left-right symmetric and hence the spatial dynamical system is reversible, but 
there is no additional requirement either for the spatial dynamical system to be conservative or for 
the temporal dynamics to be variational. In fact, many systems without these properties exhibit 
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localised states which undergo homoclinic snaking, such as the complex Ginzburg-Landau equa- 
tion with either 2:1 or 1:1 resonant forcing [22]. Similar behaviour is also observed in systems in 
higher spatial dimension, such as binary fluid convection in 2D [2J and plane Couette flow in 3D |24j . 
Therefore it is clearly of substantial interest to understand the dynamics in cases where reversibility 
is preserved but the variational/conservative structure is lost. In this paper our specific motivation 
is from work by Kozyreff and Tlidi [21] who argue that in a particular double limit of small sub- 
criticality and small critical wavenumber for the pattern forming instability, an extended version 
of the Swift-Hohenberg equation is necessary to capture the generic pattern-forming behaviour: 

dtu = rn — (1 + dxxf u + bu^ — + a{dxuf' + f3udxxU ■ (7) 

This equation is an extension of ([T|) in which, due to the double limit, quadratic nonlinearities 
containing two spatial derivatives can be established to be of the same asymptotic order as the 
other usual terms. Like ([T|), equation (|7]) is reversible; it is also equivariant under the following 
inversion involving the parameters: (x, n; r, b, a, /3) — )■ {x, —u; r, —b, —a, — /3). Equation ([7]) remains 
variational and spatially conservative for a = (3/2, but in the generic case a ^ (3/2 it loses these 
properties. On a technical level, we note that ([7|) is a semilinear PDE; extending the right-hand 
side further to include generic nonlinear terms containing four spatial derivatives would yield a 
PDE that was quasilinear but no longer semilinear. 

The dynamics and structure of localised states in ([7|) are formidably complicated. In this 
paper, we provide an initial investigation that focuses on a few key issues. Firstly, in ^we present 
the normal form analysis of the extended Swift-Hohenberg equation ^ near the pattern forming 
instability. A lengthy multiple-scales calculation allows us to compute the normal form coefficients 
explicitly in terms of the parameters in ([7]), which in turn enables us to establish the different 
regimes for the dynamics of this equation. Sections [3] and H] contain numerical investigations of 
large amplitude localised states far from onset. In ^we show that localised states in the extended 
Swift-Hohenberg equation ([7]) exhibit homoclinic snaking, though the details are not exactly the 
same as in the usual Swift-Hohenberg equation ([TJ. In 21 we examine the existence and stability 
of one-peak localised states. We find that there is both travelling wave ('drift') instability and 
standing wave instability of the single-peak state, and we discuss how these new bifurcations arise. 
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Section [5] concludes. 

2 Normal form coefficients 

In this section we use normal form analysis to examine the initial pattern forming instability of the 
extended Swift-Hohenberg equation ([7]) at r = 0. The time- independent version of this equation 
forms a fourth-order spatial dynamical system. The four spatial eigenvalues of the fixed point 
associated with the trivial state n = are given by {±(V^ - 1)1/^ ±(-^ - 1)^/2}. In r < the 
four spatial eigenvalues form a complex quartet and the origin is hyperbolic. At r = the spatial 
eigenvalues collide pairwise on the imaginary axis, and in < r < 1 they form two imaginary pairs 
so that the origin is a center. Figure [5] shows the behaviour in the spatial eigenvalues near r = 0, 
which is characteristic of the Hamiltonian-Hopf bifurcation. The dynamics in the normal form of 
this bifurcation are well understood, and the goal of this section is to classify the dynamics of the 
extended Swift-Hohenberg equation d?]) by deriving the relationship between the coefficients in the 
normal form and the parameters in ([7|). A straightforward but lengthy method to determine this 
relationship is to explicitly transform the spatial dynamical system associated with ([7]) into normal 
form. This involves introducing a four-dimensional coordinate system which reproduces the spatial 
dynamics of ([7]), then performing an appropriate linear transformation followed by a sequence of 
nonlinear near-identity transformations to match the normal form order by order. We choose an 
alternative method, based on Ref. [13] and later expanded on in Ref. [3], which involves reducing ([7]) 
to an amplitude equation and comparing this to a suitable scaled reduction of the normal form. 
The values of the normal form coefficients in terms of the parameters from ([7]) can be read off by 
comparing the two reduced equations. 

We begin in ^2.11 with a summary of the normal form for the Hamiltonian-Hopf bifurcation, 
paying particular attention to heteroclinic and homoclinic orbits because these correspond to fronts 
and localised states in ([7|) . In ^2.21 we describe the scalings that reduce the normal form ([8]) to a 
Ginzburg-Landau equation. In ^2.31 we describe a similar reduction for ([7]), and use this to derive 
the relationship between the normal form coefficients and the parameters in ([7|). In ^2.4l we describe 
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(a) r < 



(b) r = 



(c) < r < 1 



Figure 2: Spatial eigenvalues of the trivial state u = of the Swift-Hohenbcrg equation ([T]), plotted in the 
complex plane, (a) For r < the eigenvalues form two complex conjugate non-real pairs with real parts of 
equal magnitudes and opposite signs, (b) For r = the eigenvalues are purely imaginary and of multiplicity 
two. (c) For < r < 1 the eigenvalues form two purely imaginary pairs. The same eigenvalue structure 
applies to the normal form of the Hamiltonian-Hopf bifurcation (jS]) in: (a) yu < 0, (b) yU = 0, and (c) ^ > 0. 

geometrical features of this relationship. 

2.1 The Hamiltonian-Hopf normal form 

The normal form for the Hamiltonian-Hopf bifurcation is |18l [26] : 

A' = iA + B + iAP{fi;y,w) , (8a) 
B' = iB + iBP{n;y,w) + AQ{n;y,w), (8b) 

where A{x) and B{x) are complex variables, y = \A\'^ and w = -^{AB — AB), the overbar refers 
to complex conjugation, and prime denotes differentiation with respect to x in the context of our 
spatial dynamical system. The parameter /i is an unfolding parameter analogous to r in ([7]). The 
normal form is equivariant with respect to the reversibility transformation (x, A, B) — ?• (— x. A, —B). 
The terms P{fi;y,w) and Q{fi;y,w) are polynomials with real coefficients. We consider only the 
first few terms in P and Q, linear in fj, and up to quadratic order in y and w: 

P{H\ y, w) = pifi + p2y + P3W + pAy"^ + p^wy + pqw^ , (9a) 
<3(/^; w) = -qin + q2y + qsw + q^y"^ + q^wy + qew"^ . (9b) 
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Without loss of generality we can choose the sign of /i to ensure qi > 0. With this choice, the 
fixed point A = B = is hyperbolic in ^ < and a center in ^ > 0, so the bifurcation is oriented 
so that the spatial eigenvalues follow those shown in figure [2j The nonlinear behaviour of the 
unfolding near fi = depends crucially on two of the normal form coefficients, q2 and q^. The sign 
of q2 determines the criticality of the bifurcation of small amplitude periodic orbits at /i = 0: this 
bifurcation is supercritical for q2 > and subcritical for q2 < 0. The role of q^ is more subtle, as 
discussed below. We remark that if q2 is 0{1) then the role of 54 is formally irrelevant to describing 
the dynamics near fi = 0. However, the calculation that follows is performed in the neighborhood 
of the codimension-2 point (/i, 92) = (0,0) and so the sign of 94 becomes highly relevant. 

Space does not permit a complete discussion and analysis of the normal form dynamics. We 
present a very brief summary that is sufficient to highlight the differences in terms of possible orbits 
homoclinic to zero. We refer the interested reader to the analysis contained in Refs. \12 \ \17 \ \T8 \ 126) 
for further details. 

The analysis of the dynamics of the normal form is made considerable easier by the existence 
of two conserved quantities: w (defined above) and h = \B\'^ — Q{fj,; s,w)ds. Within the level 
set w = h = 0, which includes the fixed point A = B = 0, the dynamics reduces to a second-order 
nonlinear oscillator which can be conveniently written in the form: 



and by definition only y > contains physically relevant solution trajectories. 

Figure [3] indicates the shape of the potential function /(y) for the different sign combinations 
of fi, q2 and 54. In the case q^ > 0, figure [3|^a) shows that orbits homoclinic to the origin exist only 
in 52 < and only for ^um < /U < 0, where //m is defined by the condition that the discriminant of 
f{y)/y^ vanishes. At ^um the potential /(y) has a doubly-degenerate zero at a nontrivial value of y 
corresponding to a periodic orbit within the w = h = level set; there is also a heteroclinic orbit 
(i.e., a spatial front) connecting the origin to this periodic state. The codimension-one point hm is 




(10) 



where 



/(y) = Aqifiy^ - 2q2y^ - -qAy^ , 



(11) 
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(a) (b) 

Figure 3: Schematic summary of the regimes of normal form dynamics in the (/i, (72) plane, in the cases (a) 
g4 > and (b) (74 < 0, after figm-es 2 and 3 of Ref. [2^. Orbits homoclinic to the origin exist in the shaded 
regions. The dashed curves indicate A = 0, where the quantity A = (7|/4 + AqiQin/S is the discriminant of 
f{y)/y^- The sohd line in (a) labelled /zm indicates the Maxwell point. 

analogous to the Maxwell point in ([1]). While homoclinic orbits are certainly found throughout the 
shaded region in figure [3||a) , the multiplicity of localised states associated with homoclinic snaking 
is associated with the heteroclinic orbits along /iM- 

Turning to the case < 0, illustrated in figure El^b), we see that homoclinic orbits exist 
throughout fi < for q2 of either sign. The homoclinic orbits for q2 < are born at small amplitude 
as /.i decreases through zero, whereas the homoclinic orbits that exist near /i = when ^2 > exist 
at finite amplitude. Moreover, exactly at /u = these homoclinic orbits have algebraically decaying 
tails since at /U = the origin is non-hyperbolic. In addition, Dias and looss [12] have shown that, 
for §4 < and in the regime where > and q2 > 0, there exist orbits that are homoclinic to the 
periodic solutions. 
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e'B' = ie^B [e^P2\Af] + eA -e%fi + e^gsl^P + e%'- {AB - AB) + e%\A\ 



2.2 Scaling the normal form 

We wish to reduce the normal form ([8|) to an amphtude equation vahd near the bifurcation at 
H = 0. The calculations that follow are simplified by focussing on the regime of mild criticality, 
hence we restrict ourselves to the behaviour near the codimension-two point {fJ,,q2) = (0,0). 

We begin by introducing a small parameter e <C 1 and rescaling the parameters as fi = e^fi and 
q2 = e^^2- Next, we define the large spatial scale X = e^x and write {A,B) = (eyl(X), e^i?(X))e*^. 
Dropping the tildes, the polynomials P and Q in @ become 

P = eVl^l' + 0(64), (12a) 

Q = -e\i(i + e^2\A\'' + e\:,'-{AB - AB) + eVl^l' + 0{e^) , (12b) 
and so the normal form ([S]) becomes 

e^A' = e^B + ieA [e^pal^P] + ^(e^) , (13a) 

+ 0(e^). (13b) 

Z J 

Rearranging (I13ap enables us to write B in terms of A: 

B = A' -ip2A\A\^ + 0{e^). (14) 
After differentiating with respect to X, this gives: 

B' = A" - ip2 {2A'\A\'^ + A'^A') + 0(e2) . (15) 
Note that (fl3bl) gives an alternate (and independent) asymptotic expression for B': 
B' = ip2B\A\^ - qifiA + q2A\A\^ + q3^A{AB - AB) + qiA\A\^ + 0{e^) 

= ip2{A' - ip2A\A\^) \A\^ - qifiA + q2A\A\^ (16) 
+ q3'-A\A' + ip2A\A\'') - q3'-\A\\j^ - ip2A\A\'') + 94^1^1^ + 0{^) , 

after eliminating factors of B using (|14|) . Equating (|15p and (|16|) leads to a time-independent 
Ginzburg-Landau-type equation for A{X): 

A" = -qifiA + q2A\A\^ +i{p2 + ^^s) A'^' + i{^P2 - ^Q3)^Vl' (17) 
+ {q^-q3P2+pl)A\A\'' + 0{e^). 
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This equation is an approximation of the normal form ([8|) vahd in the neighborhood of the 
codimension-2 point (/i, (J'2) = (0,0). 

2.3 Scaling of the Swift-Hohenberg equation, and matching 

We now return to the extended Swift-Hohenberg equation d?]), and introduce a scahng that leads to 
an analogous version of (|17p valid near the pattern forming instability at r = 0. Again, we focus on a 
neighborhood of the codimension-2 point where the pattern forming bifurcation changes criticality. 
We anticipate that the condition of mild criticality will demand that a particular combination of 
the quadratic coefficients (6, a, /3) is small, and we allow this combination to emerge naturally in 
the calculation. 

The reduction procedure for ^ involves taking appropriate scalings for the parameters and 
expanding u{x,t) as a sum of Fourier modes multiplied by amplitudes that depend on long spatial 
and temporal scales. We introduce a small parameter e ^ 1 and rescale the parameters as: 

r = e^fi, b = bo + e^b, a = ao + e^a, /3 = /3o , (18) 

where 60, ao, /3o ^-re the values (to be determined) of the quadratic coefficients that correspond to 
q2 = 0, and fi, b, a, f3 are all 0(1). Next, we define the large spatial scale X = e^x and long time 
scale T = eH, and propose the following ansatz for solutions to ([7]): 

u(x, t) = e^G + [e^e" + e^Be^"' + e^Ce^^*^ + e^De^" + c.c] + O (e^) , (19) 

where the amplitudes 0, A, B, C, D are functions of X and T and are all 0(1), the higher 
order terms in e take the form g'^e™^ + c.c. for n > 5, and c.c. denotes complex conjugation of 
the terms preceeding it within the brackets. Note that the usual approach is to proceed order by 
order in e, bringing in appropriate modes at each order to avoid secular terms. We choose instead 
to include in (jl9p all the Fourier modes necessary, with scalings motivated in obvious ways from 
the number of quadratic interactions necessary to produce a term with the appropriate Fourier 
dependence. We proceed by substituting (jlSp and ()19p into d?]) and collecting terms with the same 
Fourier dependence e"*^, keeping careful track throughout this procedure of the size of the largest 
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neglected terms. The results for n = 0, 1, 2, 3 are: 

n = 0: O = -e2e + 5o(2eVP + e^0^ + 2e^l^n (20a) 

- 3{2e^@\A\^ + e^A^B + e^A^B) 

+ ao(2eVl^ + 8e^\Bf + 2e^{ie^AdxA - ie^{dxA)A)) 
+ /3o( - 2e^\A\^ - 8e^|Bp + 2e^(i€^{dxA)A - ie^AOxA)) 
+ 2e^h\A\^ + 2e^a|A|2 - 2e^/3|A|2 + ©(e^) 
n = 1 : e^^T^ = e^M + (^^dxxA + 6o (2e3e^ + 2e^AB + 2e^5C) (20b) 

- 3(e^A|^|2 + e^e^A + e^l^C + 2e'QAB + 2e^^|Sp) 

+ ao{Ae^AB + 12e^S(:7 + 2^^(9x9)^ + 4.ie>{dxA)B - 2ie^AdxB) 
+ /3o( - e^G^ - Be^AB - Ue^BC + 2ie^GdxA + 4ie^A5xS - 2ie^{dxA)B) 
+ 2e^6(GA + i5) + 4e^aAB - e^^{QA + BAb) + O(e^) 
n = 2 : = -Qe^S + 24ie^9x5 + 6o (e^^^ + 2e^eS + 2e^AC) (20c) 

_ 3(^40^2 ^2eVP^) 
+ ao( - e^^^ + Ge^AC + 2ie^A9xA) 
+ /3o ( - e^^^ - 4e^eS - lOe^lC + 2ie^AdxA) 
+ e4^A2 - e^dA^ - 2e^pA^ + ©(e^) 
n = 3: = -64e^C + bo{2e^AB)-{e^A^)-ao{4e^AB)-(3o{5e^AB)+0{e^). (20d) 

The n = 1 terms in ()20bp are the equivalent of the 'solvability conditions' found using other 
methods, and these are clearly the terms in which we are ultimately interested. However, we need 
first to solve (l20al) . (pOc]) and (l20dl) to express the amplitudes 0, B and C in terms of the principal 
amplitude A. We do this iteratively by writing 

G = eo + e'Ga + O(e^) , B = Bq + e^B^ + 0{e^) , C = Cq + 0{e^) , (21) 

and substituting these into (j20ap . (j20cp and ()20dp . The leading order relations are: 

Go = ci|A|2, Bo = C2A\ Co = C3^^ (22) 
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where we define the parameter combinations 



ci = 2 (60 + ao - /5o) , (23a) 

C2 = ^ {bo - ao- M , (23b) 

C3 = -J- [(26o - 4qo - 5/3o) C2 - 1] . (23c) 
64 



At next order, (|20ap and (|20cp give 



62 = ci\A\^ + C4|A|4 + ic5 {AdxA - {dxA)A) , (24a) 
B2 = C2A^ + ce^VP + iC7^5x>i, (24b) 



where we define the further parameter combinations 



C4 = 6o(c? + 2cl) - 6(ci + C2) + 8c^(ao - /3o) , (25c) 

C5 = 2(ao - /3o) , (25d) 

C6 = J [2ciC2(6o - 2/3o) + 2c3(6 + 3ao - 5^o) - 3(ci + 2c2)] , (25e) 

C7 = ^ [24c2 + (ao + /3o)] , (25f ) 



and we use the 'hat' notation in (|24|) to indicate that ci and 62 can be interpreted as higher order 
contributions to coefficients of terms that already appeared at leading order in (j22p . 

We now turn to consideration of ()20bp . As it stands, this equation contains both 0{e^) and 
0{e') terms. The problem of mixed asymptotic orders disappears if we demand that the 0{e^) 
terms vanish, and this defines the desired relation between the quadratic coefficients of d?]) at the 
codimension-2 point. Inserting ()2ip and ()22p into ()20bp . we see that the 0{€^) terms are: 

= 60 (260 A + 2ABo) - 3^1 A| 2 + AaoABo - /3o (60^ + ^ABo) , 
= [3 - ci (26o - Po) - C2 (26o + 4ao - 5/3o)] A\A\^ . (26) 
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The bracketed quantity determines the criticahty of the bifurcation at r = 0. Based on the preceding 
discussion of the normal form, we set 



q2{b,a,l3) = ^ 



3 - 2 (6 + Q - /3) (26 - /3) - - (6 - Q - /3) (26 + 4a - 5/3) 



(27) 



so that at the codimension-2 point we have (?2(^0) «0) /^o) = and the 0{e^) terms in ()20bp van- 
ish. The remaining terms in ()20bp are 0{e^) and, after substituting and tidying up, these give a 
differential equation for the principal amplitude A{X,T): 



OtA = fiA + AdxxA - 4g2^|^P + ics\A\'^dxA + ic^A^dxA + cwA\A\' 



(28) 



where: 



92 



db da 8(3 



{b,a,l3)={bo,ao,M 

cs = 6o(-2c5 + 2c7) + ao(4c7 + 2ci - 4c2) + /3o(c5 - 5c7 + 2ci + 8C2) 
eg = 6o(2c5) + ao(2ci + 4C2) + /3o(-C5 - 2C2) , 
cio = 6o(2c4 + 2c6 + 2C2C3) - 3(c? + C3 + 2ciC2 + 2c|) 
+ ao(4c6 + I2C2C3) + /3o(-C4 - 5c6 - I3C2C3) . 



(29a) 

(29b) 
(29c) 

(29d) 

Equation (]28p is the Ginzburg-Landau approximation to the extended Swift-Hohenberg equa- 
tion ([7]) valid near onset in the regimes of small criticahty — i.e., near the co dimension- two point 
(r, 92 (6, a,/?)) = (0,0). 

Both (|17p and the time-independent version of (|28p describe the spatial evolution of the ampli- 
tude A, so we can compare these two equations term by term in order to identify the normal form 
coefficients in terms of the parameters of ([7]) . This procedure gives the following: 

(30a) 
(30b) 
(30c) 

\cio ■ (30d) 



qi = 


-1/4 






P2 = 








Q3 = 








qA = 


256^ 
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For our purposes, it is sufficient to calculate only these normal form coefficients. The remaining 
normal form coefficients, as well as those associated with higher order terms that were omitted 
from Q, may be obtained by carrying out these expansions to higher order in e. 

In the limit a = (3 = 0, the above expressions reduce to the familiar normal form results for the 
usual Swift-Hohenberg equation ([I]). In this limit, becomes ^2(61, 0,0) = (27- 38b^)/36. The 
condition q2 = gives = 27/38, and at this value of b equation ()30dp gives q4\q2=o = 2202/361. 

2.4 Geometry of the parameter dependence of 52 and ^4 

In this section we describe briefly some geometrical features of the dependence of the normal form 
coefficients (72 and 54 on the parameters b, a and /3. 

Expression (|27|) for q2 is a homogeneous quadratic function of 61 , a and /3 apart from a constant 
term. After some simplification ()27p can be written in terms of a quadratic form as: 

g2(&,a,/3) = ^(27-p^Afp) , (31) 

where p = [b a (3]"^ denotes the vector of coefficients of the quadratic terms in d?]) and 

38 17 -61/2 
M= 17 -4 -17/2 (32) 
-61/2 -17/2 23 

is a 3 X 3 symmetric matrix with eigenvalues Ai ~ 66.82, A2 — 0.5113, A3 ~ —10.33. The condition 
^2 = gives 

27 = p^Mp , (33) 

which can be interpreted as an equation for a two-dimensional surface in the parameter space of 
the quadratic coefficients (6, a, /3). The signs of the eigenvalues of M indicate that this surface is a 
'hyperboloid of 1-sheet'. Figure H] shows two views of the surface (j33p . The normal form coefficient 
q2 is positive inside the hyperboloid (i.e., in the connected component of that contains the 
origin) and negative outside. The parameter inversion symmetry of d?]) is apparent in the reflection 
symmetry of the surface. 
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Figure 4: Two views of the surface 92 = plotted in the R'^ parameter space of the quadratic co- 
efficients with contours showing level sets of (74. The surface is coloured according the value 
of (74: red and blue respectively correspond to positive and negative values of (74. The open circle 
(o) marks the point {b,a,/3) = (1.8,0,0). The surface can be viewed in the accompanying movie file 
(http : //math . bu . edu/people/ jb/Research/ extendedSH/Figure_4_inovie . mpg) . 
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Figure 5: Two planar sections of the (72 = surface in figure |4j (a) a = /3/2, where equation ^ is 
variational in t and conservative in x, and (b) /3 = 0, which is the focus of ^J3] and 21 The colouring of 
the (72 = curve indicates the sign of qn, with (74 > red and (74 < blue. Solid dots (•) mark the points 
with 54 = 0. The dashed line in each frame indicates h — 1.8, and the open circle (o) marks the point 
(6, a, p) = (1.8, 0, 0). The inset in (a) shows a small ^4 > range along the curve. 

The quantity 54 is important only when 52 is small, so it is natural to consider the sign of 94 
on the surface ^2 = 0. In figure U] the regions in which 54, as given by ()30dp . is positive or negative 
are coloured red and blue, respectively. These are labelled directly in figure which also shows, 
indicated by the red and white dot, the location of the usual Swift-Hohenberg equation ([T|) for 
which a = /? = 0. This is confirmation of the result mentioned above that g4|g2=o > in ([1]). 
Moreover, it is clear from figure U] that there are substantial regions of parameter space in which 
54 < and so there is ample motivation to study the dynamics in both cases of the sign of 54. 

Recall that the extended Swift-Hohenberg equation ([7]) is variational in t and conservative in 
X when a = (3/2. This corresponds to a planar section of the three-dimensional parameter space 
defined above. Figure [5ja) shows the corresponding section of the (?2 = surface from figured] using 
coordinates (a, 5). Figure [5)||b) shows a second planar section of the (?2 = surface from figure SJ 
at /3 = 0. The behaviour of ([7]) at /? = is considered in more detail in the following section. 
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3 Homoclinic snaking in the extended Swift-Hohenberg equation 



In this section we examine the behaviour of locahsed states in the extended Swift-Hohenberg 
equation ([7]). The primary goal is simply to establish (numerically) that the homoclinic snaking 
behaviour present in ([T]) persists with the inclusion of the non-variational and non-conservative 
terms. We also point out several ways in which the localised states in ^ differ from those in ([1]). 
For simplicity, we fix 6 = 1.8 throughout this section, and use the behaviour of ([1]) at this value 
of b (as shown in figure [T]) as a point of reference. We make use of (|27p to focus on values of 
the quadratic parameters (6, a,/3) for which q2 < and \q2\ is 0{1) — i.e., the highly subcritical 
regime, where homoclinic snaking is prominent — and we use the software package AUTO [13] to 
trace out branches of localised states in the primary bifurcation parameter r. 

To begin, we examine the effect of the new terms in ([7D by increasing a from a = while leaving 
(3 = fixed. Frames (a) and (b) of figures [6] show the bifurcation diagrams of localised states at 
a = 0.1 and a = 0.5, respectively. The behaviour of stationary, even-symmetric localised states is 
qualitatively the same as that shown in figure [T] at q = 0. These localised states are organised in 
a pair of intertwined snaking branches, which we continue to label LO and LI. The saddle-node 
bifurcations on the snaking branches line up asymptotically (except perhaps the lowest two on the 
LO branch) to two r values which define the snaking region. The shaded region in each frame 
indicates the snaking region at a = 0, so increasing a increases the width of the snaking region 
and also shifts it to more negative values of r. Note however that in the case of equation ([7]) we 
can not define a Maxwell point within this snaking region. Moreover, unlike in figure [TJ we do not 
indicate the stability of any localised states in figure [6l A complete description of the stability of 
these states is beyond the scope of this article, but we do present a limited discussion of stability 
in 21 where we focus on one-peak solutions between the first and second saddle-node bifurcations 
on the LO branch. 

At the fixed values of b and /3 used in figure [6l equation ()27p for q2 has a zero at a ~ —1.306, 
with q2 < for values of a larger than this — see figure [5{b) . From ()30dp we find that > 
at this point. The normal form therefore predicts that the behaviour in the neighborhood of 
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Figure 6: Bifurcation diagrams of localised states at (a) a = 0.1 and (b) a = 0.5, with other parameters 
fixed at 6 = 1.8 and /3 = 0. The two snaking branches LO and LI contain even-symmetric, steady, spatially 
localised profiles. The profiles on the rung branches are asymmetric and drift in x at constant velocities 
±c. Temporal stability of solutions is not indicated. The shading indicates the snaking region at a = (see 
figure [IJ . The lower frames show the drift velocity of the asymmetric localised states from the lowest six 
rung branches. 
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the codimension-two point (r, a) = (0,-1.306) follows figure [3K a), just like the familiar case of 
equation ([1]) in the neighborhood of (r, 6) = (0, a/27/38). For values of a slightly greater than 
a ~ —1.306, the snaking region is exponentially narrow and close to r = 0. At larger a, the 
quantity q2 decreases to an 0(1) negative value, and the system becomes more subcritical. This 
causes the width of the snaking region to increase and shift to more negative r values. This is 
consistent with the snaking behaviour shown at the various a values in figure [H 

We note that, at the fixed values bi = 1.8 and /? = used in figure [H equation (j27p for q2 also has 
a second zero at a ~ 18.4 — see figure [5{b). The pattern forming instability at r = is subcritical 
in a G [—1.306,18.4]. Moreover, at this second codimension-two point (r, a) ~ (0,18.4) we again 
have ^4 > and so the normal form analysis suggests that nearby the behaviour of the system is 
similar to that shown in figure [3|^a) . Understanding the behaviour of ([7|) in the neighborhood of 
this second codimension-two point, and its influence on, and relationship to, the localised states 
associated with the first codimension-two point at (r, a) ~ (0, —1.306), remains an open problem. 

The bifurcation diagrams in frames (a) and (b) of figured also include the rung branches which 
cross-link the two snaking branches. The asymmetric localised states on the rungs drift in x since 
this is generic in non- variational systems. The rung branches shown in figure [6] consist of profiles 
of fixed shape that drift at constant velocity — i.e., they are travelling wave solutions of the form 
u{x, t) = u{x — ct). The drift velocity c for the various rungs is plotted in the lower frames in figure[6j 
The drift velocity varies with r along each rung, but approaches c = at the endpoints that mark 
the secondary pitchfork bifurcations where the rungs connect to the snaking branches. Recall that 
each point on a rung actually includes two different profiles related by reflection symmetry. These 
proflles drift with equal speed in opposite directions, so flgure[6Kc) includes two segments (one in 
c > and one in c < 0) for each rung in flgureEJa), and likewise for flguresEJd) andO^b). There 
are two trends in c that are apparent in figure [6l First, at fixed a the drift velocity varies from 
one rung to the next: the narrow localised states on the lower rungs tend to travel faster than the 
wider localised states on the upper rungs. Second, by comparing frames (c) and (d) it appears that 
the drift speed tends to increase with a, at least for small a. Intuitively this is to be expected as 
one moves away from the variational case q = 0. Similar trends have been observed in drift speeds 
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Figure 7: Wavcnumber k of the pattern within the locahsed states, measured along the snaking branches for 
several values of the quadratic coefficients: (i) {b, a, /3) = (1.8, 0, 0), as in figure[Tl (ii) (5, a, /?) = (1.8, 0.1, 0), 
as in figure Ela); (iii) (6,a,/3) = (1.8,0.5,0), as in figure Eljb); and (iv) {b,a,/3) = (1.8,0.1,0.2). For the 
parameter values in (i) and (iv), the system is conservative in x because a = /3/2, and the wavcnumber 
selection fc(r) is a single- valued function of r. For the parameter values in (ii) and (iv), the system is non- 
conservative in X and the wavcnumber selection k(r) is multivalued, tracing out an isola. The wavenumbcrs 
shown in the figure only apply to the wide localised states taken from far up the respective snaking branches. 

of asymmetric localised profiles in otlier systems |14j . 

The new terms in ([7]) also cause a qualitative change in the wavcnumber selection relative to 
that found in ([T|) . Figure [7] shows the measured wavcnumber k of the pattern within the localised 
states as a function of r for one complete back-and-forth cycle across the pinning region, for several 
(q,/3) pairs at fixed b = 1.8. In each case, the measurement involves wide profiles from far up 
the snaking branches where the wavcnumber is well defined. Though the measurements are made 
using one particular back-and-forth pair of segments from one particular snaking branch, they are 
characteristic of all the wide localised states far up both the LO and LI snaking branches. For the 
narrow localised states lower on the snaking branches, the wavcnumber is poorly defined. 

The curve labelled (i) in figure ([7]) shows the wavcnumber variation of the localised states along 
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the snaking branches from figure [H at a = P = 0. Recah that in this case the spatial dynamics is 
conservative so the patterns within the locahsed states must satisfy H = 0. This constraint H = 
defines the branch P of patterns included in figure [H and determines the variation in wavenumber 
k along this branch. Thus k{r) for the patterns is defined uniquely over a range of r which includes 
the snaking region of localised states; the measured wavenumber on the curve (i) in figure [7| is the 
segment of this broader k{r) curve that lies within the snaking region. In particular, we note that 
the wavenumber variation of the stable localised states along the segments of the snaking branches 
that slant 'up and to the right' is identical to the wavenumber variation of the unstable localised 
states along the segments that slant 'up and to the left'. Two profiles at the same r value on 
consecutive segments include identical patterns in their interior, and differ only in the shape of 
the fronts that connect the pattern to the flat background. Two distinct fronts are created in a 
saddle-node bifurcation at the left edge of the pinning region, and merge in a second saddle-node 
bifurcation at the right edge The H = constraint forces the fronts to approach the same 
spatially periodic orbit. 

In the absence of a spatially conserved quantity H, the saddle-node bifurcation of fronts creates 
a pair of fronts that approach two different periodic orbits. The wavenumber variation along the 
segments of the snaking branches that slant 'up and to the right' is therefore generically different 
from that along the segments that slant 'up and to the left', though of course they match at the 
saddle-node bifurcations at the edge of the snaking region. The curves (ii) and (iii) in figure [7| 
show the wavenumber measured along the snaking branches from frames (a) and (b) of figure [6j 
In these examples, the spatial dynamics is non-conservative in x and the wavenumber along the 
snaking branches traces out a loop. The larger wavenumber corresponds to segments of the snaking 
branches that slant 'up and to the right', and the lower wavenumber to segments that slant 'up 
and to the left', at least for these two examples. We note that the splitting of k{r) into a loop has 
also been observed in other spatially non-conservative systems, such as plane Couette flow |14j . 

The curve labelled (iv) in figure [7] shows the wavenumber variation at a = 0.1 and /? = 0.2. 
As a = (3/2, this corresponds to a spatially conservative limit of d?]). As such, the wavenumber 
variation is again a single- valued function of r. The splitting of the wavenumber selected by the 
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localised states is therefore a measure of how spatially non-conservative is the system, though the 
high-precision numerical measurements required to observe the splitting may make this impractical 
as a diagnostic tool. 

4 One peak solutions 

In this section we investigate another aspect of how the non- variational and non-conservative terms 
effect the behaviour of d?]). We focus on the stability of one-peak localised states since this enables 
us to show clearly the emergence of two new kinds of temporal instabilities: drift and standing 
oscillations. The one-peak localised states occur on the segment of the LO branch between the first 
and second saddle-node bifurcations. It is convenient to introduce the notation LOipk to refer to 
this segment. We organise the results by fixing b = 1.8 and /3 = 0, and examining the states on 
LOipk at increasing values of a. 

Throughout this section, we report on the linear stability of the one-peak states. As usual, 
the linear stability of a stationary solution uo{x) to ^ is determined by substituting u{x,t) = 
uo{x) + eU{x)e'^'' into ([7|) and ignoring 0{e^) terms. The mode U{x) and growth rate a satisfy the 
linear eigenvalue equation aU = C[uq, dx]U where the linearised operator associated with the right 
side of d?]) is: 

C[uq, dx\=r-{l + dxxf + 26tto - Suq + 2a{dxUo)dx + /3 {{dxxUo) + uod^x) ■ (34) 

Spatial translation invariance of ([7]) implies that any steady solution always has a neutrally stable 
((T = 0) Goldstone mode Ug{x) = dxUo{x). The one-peak solutions of interest here are even- 
symmetric, and the associated Goldstone mode is odd-symmetric. In what follows, we report the 
existence of unstable {Re a > 0) modes using the notation [ ], where (fn-c) is the 

number of even-symmetric modes with real (complex) eigenvalues a, and tLj. (uc) is the number 
of odd-symmetric modes with real (complex) eigenvalues. Stable solutions are labelled [0,0,0,0]. 
Note that both nic and iic must be even; in what follows, we never observe unstable odd-symmetric 
modes with complex eigenvalues so nc = throughout. 
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Figure 8: Bifurcation diagram of one-peak states for {b,a,(3) = (1.8,0.5,0). (a) Snaking diagram of 
localised states, with no indication of stability, (b) Detail of the LOipk segment; solid (dashed) curves 
indicate stable (unstable) states. The number and symmetry of unstable modes is also indicated, using 
the notation defined in the text, (c) Detail of the right edge of LOipk- Solid dots (•) indicate saddle-node 
bifurcations; the open circle (o) indicates the secondary bifurcation to the rung branch. One-peak solutions 
are stable over the entire LOipk segment. 




Figure 9: One-peak profiles and eigenfunctions from figure H] at (6,a,/3) = (1.8,0.5,0), (a) near the left 
edge of LOipk and (b) near the right edge of LOipk. Each column includes the stationary one-peak state u(x), 
the neutrally stable Goldstonc mode Ug{x), and the even-symmetric amplitude mode Us,{x) that changes 
stability at the corresponding saddle-node bifurcation. The right column also includes the odd-symmetric 
phase mode Up{x) which is associated with the lowest rung branch. 
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A natural starting point is the case q = 0, as the behaviour of the usual Swift-Hohenberg equa- 
tion ([1]) is well-known. As shown in figure [H the one-peak localised states on LOipk are all stable. 
Increasing a initially maintains qualitatively very similar behaviour. Figure [8] shows the behaviour 
at a = 0.5. Figure [S[a) shows the usual snaking behaviour of LO and LI. Figure [5[b) shows a 
detail of the LOipk segment, and includes stability assignments. Below the saddle-node bifurcation 
at the left edge of LOipk, ^ single even-symmetric 'amplitude' mode becomes unstable. Likewise, at 
the saddle-node bifurcation at the right edge of LOipk an even-symmetric amplitude mode becomes 
unstable. This is followed by an odd-symmetric 'phase' mode which becomes unstable slightly 
above this saddle-node bifurcation, at the secondary bifurcation that gives rise to the lowest rung 
branch. Figure O shows the profiles u{x) at two points along LOipk, one near the left edge and 
one near the right edge. This figure also includes plots of the modes that are most important in 
determining stability. We note in passing that, unlike in the case of the usual Swift-Hohenberg 
equation ([T]), the amplitude mode which passes through zero growth rate at the left edge of LOipk 
is different to the one that passes through zero growth rate at the right edge. The former merges 
with another even-symmetric mode along and becomes complex with Recr < 0, corresponding to a 
stable oscillatory mode. 

Qualitatively new behaviour occurs at a = 1.6, as shown in figures [TO] and [TTJ The bifurcation 
behaviour in the neighborhood of the left edge of LOipk remains unchanged — below this saddle- 
node bifurcation there is a single unstable amplitude mode which stabilises at the saddle-node 
bifurcation. However, the profiles on LOipk are not stable all the way up to the saddle-node 
bifurcation at the right edge. Instead, a phase mode loses stability slightly before the right edge. 
This is the same instability as that which was previously associated with the secondary bifurcation 
to a rung branch, so now the rung connects to the LO branch slightly below the right saddle-node 
bifurcation. Though not labelled in the figure, we find that the profiles on the rung branch remain 
unstable. 

Increasing a causes the secondary bifurcation to the rung branch to move further from the right 
edge of LOipk, but also leads to new instability. The bifurcation diagram at a = 2.0 is shown in 
figure [T2j At this value of a, the first mode to become unstable along LOipk is an even-symmetric 
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Figure 10: Bifurcation diagram of one-peak states for (6, a,/3) = (1.8,1.6,0). Notation follows figure [H 
Stationary one-peak solutions are stable from the left edge of LOipk up to r ~ —0.2358 where the odd- 
symmetric mode associated with the rung loses stability. 




Figure 11: One-peak profiles and eigenfunctions from figure [TOl at = (1.8,1.6,0), (a) near the left 

edge of LOipk and (b) near the right edge of LOipk- Notation follows figure O 
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Figure 12: Bifurcation diagram of one-peak states for {b,a,f3) — (1.8,2.0,0). Notation follows figure [5] 
Stationary one-peak solutions are stable from the left edge of LOipk up to r ~ —0.2682 where an even- 
symmetric oscillatory mode loses stability; the oscillatory instability is marked with a diamond symbol (o). 




Figure 13: (a) One-peak profile and eigenfunctions from figure[T2]at (6, a, (3) = (1.8, 2.0, 0), at a value of r 
slightly above the oscillatory instability on LOipk. Notation follows figure El The eigenfunction Uosc{x) of the 
mode associated with the oscillatory instability is also included; solid and dased lines indicate the real and 
imaginary parts of Uosd^) respectively, (b) Space-time plot of an oscillon at (6, a, /?) ~ (1-8, 2.0, 0), shown on 
the domain (x, t) E [—50, 50] x [0, 100]. The r value for this solution is r = —0.2681, which is slightly above 
the oscillatory instability of the one-peak solutions from figure [T^l The oscillon can be viewed in the accom- 
panying movie file (http://math.bu.edu/people/jb/Research/extendedSH/Figure_13_movie.mpg). 
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oscillatory mode, which loses stability at r ~ —0.2682. The relevant mode Uosc{x) is therefore 
complex, and its real and imaginary parts are shown in figure [TSTa). This mode eventually re- 
stabilises much further along LO, above the saddle-node bifurcation at the right edge of LOipk. We 
find that the imaginary part of the complex eigenvalue at onset is Imo" ~ 0.55 and that it remains 
nearly constant along the snaking branch. 

The stationary localised solutions on LOipk at r values above the complex instability are un- 
stable, but may evolve to time-dependent states which remain spatially localised. Such solutions 
are often referred to as oscillons. Figure [TSTb) shows the space-time plot of one such solution at 
r = —0.2681, slightly above the initial instability on LOipk at r ~ —0.2682. The criticality of this 
instability is, however, difficult to confirm. Numerical results indicate that the oscillons that occur 
above onset are large perturbations of the unstable stationary solutions on LOipk. Furthermore, 
the oscillons persist as r decreases below r ~ —0.2682. One explanation for this might be that the 
oscillatory instability is subcritical, and that the branch of oscillons that emerges from this point 
turns around and stabilises in a saddle-node bifurcation. A second possibility is that the branch 
that emerges from LOipk remains unstable, and the localised solution shown in figure [TSTb) lies 
on a separate branch of localised oscillations that emerges elsewhere in the bifurcation diagram. 
This second possibility occurs in the autonomous system of reaction-diffusion equations studied in 
Ref . [1] , which includes a bifurcation structure very similar to that shown in figure [T2l 

Further increase in a causes the oscillatory instability to invade more of the stable range of the 
LOipk segment. At a ~ 2.8 this instability reaches the left edge of LOipk so all stationary one-peak 
solutions are unstable. 

5 Conclusions 

In this paper we have considered the effect of non-variational and non-conservative terms on the 
well-known quadratic-cubic Swift-Hohenberg equation. Such an investigation is well motivated 
in general by the wealth of recent work on homoclinic snaking in a variety of contexts where no 
variational principle exists, and in particular by the paper of Kozyreff and Tlidi [21] who show 
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that ([7]) arises naturally as a model equation for a long-wavelength pattern forming instability in 
one-dimensional dissipative systems. 

The dynamics of the extended Swift-Hohenberg equation ([7|) are extremely rich and we have 
been content to set the scene for future work by focussing on a few aspects of the problem. Firstly we 
have carried out the computation of the normal form coefficients for the pattern-forming bifurcation 
problem. This is a key first step in the analysis of this extended equation. It allows us to use the 
normal form results as a guide to the kinds of localised states that should exist in different parameter 
regimes. One important result is the existence of substantial open regions of parameter space in 
which q4 < 0. 

Next we investigated numerically the effect of one of the new terms on the snaking behaviour 
at large amplitude, far from the initial linear instability. We found that the snaking phenomenon 
persists, and we pointed out aspects of snaking in the extended Swift-Hohenberg equation ([7]) that 
differ from the well-known snaking in ([1]), including wavenumber selection within the localised 
states, and the behaviour of the solutions from the rungs. 

Lastly we investigated the emergence of two new kinds of instability of stationary one-peak 
localised states that affect the parameter range over which stable one-peak states exist. These 
are (i) a drift instability and (ii) an oscillatory instability. The latter may lead to stable localised 
oscillations, but the bifurcation structure of such solutions remains unclear. We showed how the 
two new instabilities arise as qualitative changes to the traditional snaking bifurcation diagram. 
Since both instabilities are generic we expect that this analysis will provide a guide to the location 
and dynamics of instabilities of multi-peak localised states and quite possibly multi-pulse localised 
states as well. The details of these cases, as so much else, we leave for future work. 
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